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Thermoacoustic engines convert heat energy into acoustic energy. Then, the acoustic energy can be used 
to pump heat or to generate electricity. It is well-known that the acoustic energy and therefore the exer¬ 
getic efficiency depend on parameters such as the stack’s hydraulic radius, the stack’s position in the 
resonator and the traveling-standing-wave ratio. In this paper, these three parameters are investigated 
in order to study and analyze the best value of the produced acoustic energy, the exergetic efficiency 
and the product of the acoustic energy by the exergetic efficiency of a thermoacoustic engine with a par¬ 
allel-plate stack. The dimensionless expressions of the thermoacoustic equations are derived and calcu¬ 
lated. Then, the Particle Swarm Optimization method (PSO) is introduced and used for the first time in the 
thermoacoustic research. The use of the PSO method and the optimization of the acoustic energy multi¬ 
plied by the exergetic efficiency are novel contributions to this domain of research. This paper discusses 
some significant conclusions which are useful for the design of new thermoacoustic engines. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

The thermoacoustic energy conversion is interpreted as the 
thermodynamic interaction of a compressible moving oscillating 
fluid with a solid wall which has a temperature gradient in the 
direction of the oscillation. During the displacement of the fluid, 
the variations of the pressure within the fluid generate heat trans¬ 
fers with the wall of the solid, constituting the stack or the regen¬ 
erator of the thermoacoustic device. According to whether the fluid 
is heated or cooled during the phases of compression or expansion 
and according to the direction of its displacement, the fluid carries 
out a thermodynamic cycle which makes it possible to generate or 
consume mechanical/acoustic energy. When this energy is gener¬ 
ated, the thermoacoustic device is called an engine or a prime 
mover, and when the energy is consumed, the thermoacoustic de¬ 
vice is called a refrigerator or a heat pump. In the case of the ther¬ 
moacoustic engines, the energy generated can thereafter be used to 
drive loads such as electrical generators or pulse tube refrigerators. 

In fact, the thermoacoustic devices have many advantages. 
Firstly, these devices can use any external energy sources like solar 
energy, radio isotope sources, combustion of any fuel oil, biomass, 
hot gas downstream from cycle, etc. Secondly, these devices do not 
contain any moving part except the driven loads in the devices 
coupled to the thermoacoustic systems (for example, the electric 
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production devices or the actuator of pressure used in the cooler 
devices). Flence, it is obvious that the thermoacoustic systems offer 
a high level of reliability and low maintenance requirements. Fur¬ 
thermore, as shown in Fig. 1, these devices basically consist of a 
stack (when using a standing-wave) and/or a regenerator (when 
using a traveling-wave), a cold and a hot heat exchanger, and a res¬ 
onator. That is why the realization cost of these thermoacoustic 
systems is quite low comparing to other conventional systems. 
And last, but not least, these devices are also friendly environmen¬ 
tal because they use a non-toxic working gas (air, helium, nitrogen, 
etc.). 

The interest to thermoacoustic systems has been arisen and 
expanded since the early 1980s when Swift from Los Alamos 
National Laboratory and Garrett from PEN State University, whom 
can be regarded as the precursors in the field [1,2], have success¬ 
fully worked and made the first practical thermoacoustic devices 
that produced a useful work. 

Based on the works of Rott [3-9], Swift [10-12] has shown that 
the boundary functioning between thermoacoustic stack-based en¬ 
gines and heat pumps corresponds to a critical mean temperature 
gradient across the stack. When the value of the wall temperature 
gradient is higher than the value of the critical mean temperature 
gradient, the acoustic wave is created and the thermoacoustic de¬ 
vice works as an engine (a prime mover). Moreover, when the va¬ 
lue of the wall temperature gradient is less than the value of the 
critical mean temperature gradient, the acoustic work is consumed 
and the thermoacoustic device functions as a refrigerator (a heat 
pump). 
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Nomenclature 



A 

cross-sectional area (m 2 ) 

y 

ratio, isobaric to isochoric specific heat 

c 

speed of local sound (m s -1 ) 

3 

penetration depth (m) 

Co 

speed of sound (m s -1 ) 

£ 

porosity 

C P 

isobaric heat capacity (J Kg -1 K -1 ) 

Es 

plate heat capacity ratio 

f 

frequency (Hz) 

P 

density (kg m -3 ) 

H x 

total energy flux in the x direction (W) 

A 0 

wave length (m) 

ko 

wave number (m -1 ) 

Ilex 

exergetic efficiency 

k 

thermal conductivity (W m _1 I< -1 ) 

Y]en 

energetic efficiency 

keq 

equivalent thermal conductivity (W itT 1 K _1 ) 

p 

dynamic viscosity (kg m -1 s -1 ) 

L 

stack’s length (m) 

T 

traveling-standing-wave ratio 

P 

pressure (Pa) 

CO 

angular frequency (rad s -1 ) 

Pr 

Prandtl number 



Qx 

heat flux in the x direction (W) 

Subscripts 

r 

gas constant (J Kg -1 1C 1 ) 

0 

ambient 

r h 

hydraulic radius (m) 

a 

first order oscillation 

S 

entropy flux (W K -1 ) 

g 

gas 

T 

temperature (K) 

s 

solid 

t 

time (s) 

t 

thermal 

u 

velocity (ms -1 ) 

V 

viscous 

u 

volume flow rate (m 3 s -1 ) 



w x 

work flux in the x direction (W) 

Superscripts 

0 

time-averaged 

f 

mean, space-averaged perpendicular to x of/ 

»[] 

real part of 

f 

time rate of/ 

3[] 

imaginary part of 

f 

conjugate part of/ 



r 

dimensionless number of/ 

Greek letters 



oc 

thermal expansion coefficient (K ') 




In addition to the two classes of thermoacoustic devices (en¬ 
gines and refrigerators), we can also classify each thermoacoustic 
devices as a standing-wave or traveling-wave (Fig. la and b respec¬ 
tively). In a standing-wave (stack-based devices), the phase be¬ 
tween the oscillating velocity of the gas and the oscillating 


(a) 




Fig. 1 . General structure of a thermoacoustique engines: (a) standing-wave; (b) 
traveling-wave. 


pressure is close to 90°, and the gas interacts with the stack’s wall 
by a voluntarily imperfect thermal contact. While in a traveling- 
wave (regenerator-based devices), the phase is close to 0°, and 
the gas and the regenerator’s wall are in very good thermal contact. 
In 1979, the American engineer Ceperley [13] has discovered that 
the Stirling machines are not other than traveling-wave thermoa¬ 
coustic devices. In such a device, energy conversion efficiency 
can approach the theoretical maximum efficiency, i.e. Carnot effi¬ 
ciency, which depends only on the temperatures of both the hot 
and the cold thermal reservoirs. The prototypes achieved by Back- 
haus and Swift [14-18] of “thermoacoustic” Stirling machines have 
opened the way to the realization of increasingly powerful sys¬ 
tems. The mini electric thermoacoustic generator produced by 
Backhaus et al. [19,20] for NASA offered already an electric output 
power of 58 W with a total efficiency of 18%, thus approaching the 
25% of the efficiency of the gasoline engines with internal combus¬ 
tion (efficiency of the sink to the wheels). 

At the present time, thermoacoustic systems, either the engines 
or the refrigerators, do not reach sufficient exergetic efficiency (the 
ratio between the system conversion efficiency and the theoretical 
Carnot efficiency). The most values of output exergetic efficiency 
usually published are between 10% and 20%, and at best 30%, but 
an objective of 40% seems however to be realistic if the various 
sources of energy dissipations or degradations are identified and 
decreased. For example, in the pulse-tube refrigerators, one can 
avoid using devises, such as long capillaries, in order to realize 
the phasing between the pressure and volumetric velocity oscilla¬ 
tions, because these are acoustic impedances, which represent a 
degradation source for the acoustic power. However, one can re¬ 
place these impedances by adding a feedback acoustical loop to 
the acoustic network in order to recycle a large part of the acoustic 
energy carried out from the regenerator back into its entry. But, by 
doing this, one open the door to other categories of acoustic energy 
degradations which are named “streaming flows”. However, it is 
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possible to decrease this flows by introducing additional bodies 
into the acoustic network (for example, jet pump dissymmetrical 
membranes, openings, conical tubes, etc.) but these elements 
cause, also, a pressure drop and friction dissipations of acoustic 
power. 

Systematic studies on the causes and consequences of the en¬ 
ergy losses, using the influence of the traveling-standing-wave ra¬ 
tio, become necessary to realize and design high efficiency 
thermoacoustic devices. Recently, it appeared in the specialized lit¬ 
erature, some studies related to the influence of the traveling- 
standing-wave ratio or treating the energetic optimization of ther¬ 
moacoustic systems. Kang et al. [21-23] have studied the thermoa¬ 
coustic effect of traveling-standing-wave and they have done an 
optimization study of hydraulic radius of a thermoacoustic engine. 
Very often, the results of the author’s works are based on the exclu¬ 
sive framework of a pure traveling-wave or standing-wave with a 
sound field presumably given within the resonator and not modi¬ 
fied by the presence of the stack and the heat exchangers. Hu et al. 
[24] have studied the evaluation of thermal efficiency and energy 
conversion of thermoacoustic Stirling engines which is working 
with just a pure traveling-wave. Babaei and Siddiqui [25] have 
developed an optimization algorithm to design thermoacoustical- 
ly-driven thermoacoustic refrigerators using a pure standing-wave 
thermoacoustic design. Yu and Jaworski [26] have studied the im¬ 
pact of acoustic impedance on the acoustic power produced in pure 
traveling-wave of thermoacoustic devices. 

This paper tries to introduce more advances in the thermoacou¬ 
stic research by indicating how the stack’s hydraulic radius, the 
stack’s position in the resonator, and the traveling-standing-wave 
ratio can modify the acoustic energy, and therefore the exergetic 
efficiency and most importantly the product between the acoustic 
energy and the exergetic efficiency of a thermoacoustic engine 
with a parallel-plate stack. For doing that, the thermoacoustic 
equations, developed by Rott and written by Swift [11], are rewrit¬ 
ten in a dimensionless numbers in order to facilitate the calculus 
and the analyses. Next, synthetic optimization is represented to 
give a schematic description of the exergetic efficiency vs. the 
acoustic energy. Then, the Particle Swarm Optimization method 
(PSO) is introduced and used. So far in the literature, synthetic 
optimizations have been applied to the thermoacoustic research. 
Also, classical optimization methods such as the gradient-descend 
and the quasi-newton cannot be involved easily in the thermoa¬ 
coustic optimization problems, because these ones are so compli¬ 
cated, irregular and non-linear. Consequently, classical and 
synthetic optimization methods need a huge amount of calculation 
time and space in the actual computers in order to optimize a func¬ 
tion depending on more than two parameters, which can be fre¬ 
quently faced in the thermoacoustic research. In contrast, the 
PSO method is a powerful optimization method that reduces the 
calculation time and space and that could be easily used. 

As a second novel contribution, the optimization of the product 
of the acoustic energy and the exergetic efficiency is studied. This 
leads to better results because it finds the best tradeoff between 
the two parameters of the product. This paper is ended up with a 
discussion of results and a significant conclusion that could be 
used as a guide to design any new thermoacoustic devices. 

2. Acoustic power amplification and attenuation in 
thermoacoustic devices 

The acoustic wave amplification was explained qualitatively a 
long time ago by Rayleigh [27]: “if heat be given to the air at the mo¬ 
ment of greatest condensation [i.e., greatest density] or be taken from 
it at the moment of greatest rarefaction, the vibration is encouraged". 
Actually, this phenomenon occurs only when an acoustic wave 
oscillates in a fluid in contact with a porous medium (stack or 


regenerator), which has a sufficient temperature gradient. The 
time-averaged acoustic power W over a cross-sectional area “A" 
of the channel is calculated as: 

(W) = {P a ha) =^\p a U a ] ( 1 ) 

where p a is the oscillating pressure, ii a is the oscillating volume flow 
rate, which is the integral of the x components of the oscillating 
velocity u a over the cross-sectional area A of the channel, 5R[] is 
the real part of the term between the brackets, () denotes the 
time-average, the overdot denotes the time rate, the overbar de¬ 
notes the spatial-average perpendicular to x, and the tilde denotes 
the complex conjugate. 

The enthalpy flow is generally introduced to explain the physics 
and the dynamics of traditional Stirling machines, of pulse-tube 
refrigerators or of thermoacoustic devices. In fact, it is considering 
the first law of thermodynamics, which is called the conservation 
of energy principle. So, the acoustic approximation to the time- 
averaged total power (H) over a cross-sectional area A can be ex¬ 
pressed as: 

(H) = \c K p g R[T a ii a ] = {p a u a ) + T g (S) = (VV> + (Q) (2) 

where (S) is the time-averaged entropy flux over a cross-sectional 
area A, C pg is the fluid heat capacity at constant pressure, p g is the 
fluid mean density, T a is the oscillating temperature amplitude of 
the fluid, T g is the fluid mean temperature and (Q) is the time-aver¬ 
aged heat power over a cross-sectional area A. As shown by Eq. (2), 
the total power is equal to the sum of the acoustical (or mechanical) 
power and the heat power. 

From Eq. (2), we can distinguish that when the fluid is oscillat¬ 
ing within a porous media, characterized by its hydraulic radius r h , 
which is much smaller than the thermal penetration depth S t 

(where <5 t = I, which is the case of nearly ideal regenera- 

V V a '"s c p / 

tors, the oscillating temperature in the fluid is almost null T a ~ 0, 
even if the oscillating pressure is not null. This means that the total 
power (H) is also null and the acoustic power (VV) is always equal 
and opposite to the heat power (Q). In a refrigerator or a heat 
pump, the consumption of the acoustic power (W) induces the 
existence of a heat flow (Q) that goes from the cold zone towards 
the hot zone. On the other hand, when the thermal penetration 
depth is much smaller than the hydraulic radius, and if there are 
no radial losses of heat in the fluid, so the system is in adiabatic sit¬ 
uation and the total power remains constant. 

As said previously, thermoacoustic devices can be classified as 
standing-wave stack-based devices or traveling-wave (acoustic- 
Stirling or pulse-tube) regenerator-based devices. Indeed, it should 
be remembered that a pure standing wave is a combination of trav¬ 
eling-wave and reflected-wave at same frequencies. 

The regenerators and stacks must be sandwiched between two 
heat exchangers at different temperatures to work usefully. In fact, 
there is an intrinsic difference between the stack-based and regen¬ 
erator-based devices. In the first one, the stack’s hydraulic radius 
must be nearly equal to the thermal penetration depth (r h ~<5 t ), 
in order to provide an imperfect thermal contact between the fluid 
and the stack’s wall. This condition is needed to create the thermo¬ 
dynamic cycle, but the temporal phasing between the oscillating 
velocity and the oscillating pressure is close to 90°, which might 
carry no acoustic power and subsequently the efficiency is rela¬ 
tively low. In such devices, when the temperature gradient is high¬ 
er than a critical mean temperature gradient, the stack produces an 
acoustic power without needing the presence of an initial acoustic 
power to create more acoustic power. In addition, the power 
producing phenomena in the stack is not altered by the presence 
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of acoustic power flowing through it in the either direction. Other¬ 
wise, when the temperature gradient is less than the critical mean 
temperature gradient, the acoustic power, flowing in any direction 
into the stack, is absorbed. 

In the regenerator based devices, the hydraulic radius is much 
smaller than the thermal penetration depth (r h -c <5 t ) in order to 
satisfy the excellent thermal contact between the fluid and the 
regenerator’s wall. This condition makes approximately an ideal 
thermodynamic cycle, so the temporal phasing between the oscil¬ 
lating velocity and the oscillating pressure is close to 0° and subse¬ 
quently the efficiency of such devices is better than the stack- 
based devices. But the regenerator only amplifies or attenuates 
an existing acoustic power in it, so if there is no acoustic power 
incoming in the one end of the regenerator, there will not be an 
acoustic power outgoing in the other end. The acoustic power is 
amplified only if it flows from the cold to hot end of the regenera¬ 
tor, and the acoustic power amplification factor is equal to the ratio 
of the hot to cold temperature end of the regenerator. Contrariwise, 
the acoustic power is attenuated or absorbed in the regenerator, 
when it flows from the hot to cold temperature end of the regen¬ 
erator. Because the high efficiency of regenerator-based devices, 
it is interesting to realize an acoustic network with a feedback 
acoustical loop with the possibility of modifying the traveling- 
standing-wave ratio. 


3. Dimensionless of thermoacoustic equations 


Since the length of the stack is too small in comparison to the 
wave length, the basic structure of the acoustic pressure field is as¬ 
sumed to be not fundamentally modified by the presence of the 
stack. Hence, the oscillating pressure propagating, in an ideal gas 
in the direction x of a channel with a cross-sectional area A , is sup¬ 
posed to be: 

p a (x)e> mt = (1 - x )p_ + 1 p + = [p(( 1 - z)e ik ° x + ze^ k ° x )]e ialt (3) 


where r is the traveling-standing-wave ratio, 0 ^ t ^ 1, it is equal 
to 1 for a pure traveling-wave, 0.5 for a pure standing-wave, and 
0 for a pure reflecting traveling-wave, p a is the non-time dependent 
part of the oscillating pressure, p is the pressure amplitude, k 0 is the 
wave number, k 0 = g and c 0 is the sound speed at the refer¬ 
ence temperature T 0 , c 0 = \JyrT 0 . 

So, the pressure of a pure-standing wave, in an ideal acoustic 
channel, is written as: 


P a (x)e imt = ^P +\Pr = 


-p(e ik ° x + e~ J * oX ) 


(4) 


And the pressure of a pure and reflecting traveling-wave, in an ideal 
acoustic channel, are respectively expressed as: 

f PflMe"* = P+ = \p(e-i k ° x )]e> wt 
| p a (x)e'“' t — p_ = \p(e' k ° x )]e iwt 

Next, the thermoacoustic equations derived by Rott and written 
by Swift [11] are rewritten in a dimensionless number. 


3.1. Acoustic field 

The pressure field supposed in Eq. (3) is rewritten by using def¬ 
initions in Table 1, in a dimensionless number as: 

p* = DR{{ 1 - z)e> 2nx " + ze- j2nx ‘), (6) 

where DR is the drive ratio, DR = f . 

Ps 

The Rott’s acoustic approximation to the x components of the 
momentum equation gives the relation between the oscillating 


Table 1 

Dimensionless number definitions of parameters. 


Parameters 

Parameters in dimensionless number 

Position, x 

x * = J 

Length of the stack, L 

L* — — 

L 

Pressure, p a 

n* = E° 

Pg 

Velocity, u a 

U * =IHa 
u a c 0 

Temperature, T h \ T g ; T a 

T* — Th- T* _ T g' T* — 

- V l g ~ V ‘a~T^ 

Density, p g 

Ta 

II 

Work flux, (VV X > 

/w*\ = y(m. 

\ */ 75 Pg^0 

Work flux gradient, 

"W _ yflfr.) 

dx’ ApgCo 

Heat flux, (Q*) 

/n*\ = m*- 

\ V ' X / ApgCo 

Energy flux, (H x ) 

\ */ TWi 


velocity and the oscillating pressure gradient. This relation in a 
dimensionless number becomes: 


dp* a = ojA q p; 

9x* jc 0 (1 -g„(s.)) 


(7) 


where g 0 (s») is the spatial averaged of a complex thermoacoustic 
function which depends on the specific channel geometric for two 
infinite parallel plates, g 0 (s,) — tan s h(s * ) ,s, =^y/2j, where 5 V is the 
viscous penetration depth, <5 V = y/Jt, Pc is the Prandtl number, 
Pr = 4 and p g is the mean density for an ideal gas, p g = -St- — — ^ 

o t rT g Cpg(y-1)Tg 

Transforming the Rott’s acoustic approximation to the continu¬ 
ity equation into a dimensionless number, we deduce that: 


dul 
dx * 


CD/l 0 
Fo 


■ (7 - 1 )&> (s. ■ y/Pr) ] Pa + Go i ^ K 


( 8 ) 


where C 0 is a combined of the spatial averaged thermoacoustic 
functions, 

_ (go(s.)-£o(s.v / Pr)) 

Go ~ (Pr - 1)(1 -g 0 (s„)) ' 


3.2. Acoustic power 

Because we are just interested by what is happening inside the 
stack, the studying zone is reduced to the stack only, given that the 
position of the studying zone varies in the resonator. So the posi¬ 
tion in the stack, x*, varies between [x* -|;X*+f], (see Fig. 2), 
where x* is the center of the stack, and L* is the length of the stack. 

Under these circumstances, the acoustic power at position x* is 
equal to: 

(w* x ) = 2 yi [ p “ !J “] = ^[PXj (9) 



1 



| 





1 111 1111111II11 
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Fig. 2. The studying zone and the stack’s position in the resonator. 
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The derivation of the time-averaged acoustic power d(W*j pro¬ 
duced in a length dx* of the channel is given as: 

d(vv*\ i dTt / . \ l dT* r , 


^0?Pg| ,|2 W„, |2 

\ra\ n 


2 c 0 r f 


2 yp g 


( 10 ) 


where T* is the gas mean temperature which is supposed to be lin¬ 


ear throughout the stack, TJ = ^ = T h + - 


- (1 - T* h ), and its ax¬ 
ial thermal gradient ^ J h \ where r v and i are the acoustic 
viscous resistance and the acoustic thermal resistance, respectively, 

3(-go(s.)) 

1-goM 2 rp 


T v = fflftr 


and 


3.3. Heat power 


The heat flux, <^Q*y, which is a combination of two terms, an 
internal heat flux and a thermal-conduction, is equal to: 


<&) = (internal 
&) = 


Co 1 


2m/o Pg 

yT 0 

PgCa'-o 


Qcond 

d P: 


dx* 


Pa(So(S t ) - lA) 


ypJo h * 


C 0 )? Q eq PI dx* 


, ,, dT! , dT* 

(l- e),<s # + ^d# 


( 11 ) 


( 12 ) 


3.4. Total power 


where N is the discretization number of the stack, x’t is the center 
position of the ith interval step of the discretized stack, and ^ is 
expressed in Eq. (10). 

Supposing that a system receives a quantity (Ql) of heat from a 
heat source which has an imposed hot temperature T* h , and the sys¬ 
tem produces a quantity -d(VV*) of acoustic power, then the en¬ 
ergy efficiency of the system is: 



(16) 


And therefore the exergetic efficiency is defined as: 


rien -r h A(W') 

Carnot {T^ ~ 1) (0_l) 


(17) 


where i] C amot , is the Carnot efficiency, (Q0 is the heat added to the 
system and can be calculated by using Eq. (12): 



( Co 1 _ 

y2c0/lo PI J 

(18) 


dPl 

dx’ 


Pa(go(S,) - I/O 


+m°K 

Co An tH P: 


1 at: 


dx * 


where x* h is the position of the hot heat exchanger in the resonator. 

Replacing Eqs. (15) and (18) into Eq. (17), the exergetic effi¬ 
ciency becomes: 


I?ex 


-n 




(n-i) 


At* 3 SwC&M-h 


1 i ypg T p i/ j_ dr l 


(19) 


For an ideal gas and for a material with high heat capacity, the 
total power is the sum of the acoustic power and heat power 
according to the first law of thermodynamics, and is expressed as: 

(13) 


h;>=TO+(q ; 


h: )= 


Co r- 


2<DAoPg 

yT 0 


dPa 

dx' 


PaO-i) 


ypJo k < i dT g 

CqAq 


,3 >*dx* 


PgCo^o 




(14) 


where 




keq 2m 3 p 0 (l-Pr) 


= Pf-’ and i> = So(s,) + 

ri 0 


dp a 


dx* 


pg T o 1 , 

,pt= — = ^ = ^ and p 0 
Po T g Tt Fo 


So(s,VPr) -So(s.)) 
(1 +Pr) 


4. Exergetic efficiency of a thermoacoustic engine 


Approximately, the acoustic power produced or consumed by a 
stack or a regenerator can be calculated as the sum, throughout the 
stack, of the derivation of the time-averaged acoustic power 
d(VV*\ produced in a length dx*, so: 


. . ,*Nd W: 


i= 1 

1 dT: 


dx* 


1=1 \ g 


,_ Llii-r ura*r,*i ^oyPg \„*\ 2 2oC 0 r„ U|t| 

+ 2Tjdx* /[ Gol/ ^ a 2c 0 r t ‘ Pa ' 2yp g I “I 


(15) 


5. Synthetic optimization 

In this part, a synthetic optimization is represented and dis¬ 
cussed. The idea that lies behind synthetic optimization is to deter¬ 
mine how to modify the exergetic efficiency as a function of the 
acoustic energy according to different values of r h , z and x* while 
all the other parameters remains fixed and according to two differ¬ 
ent working gases Helium and Air. The mean pressure is fixed at 
25 bar, the dimensionless of the stack’s length L* is equal to 
0.0025 and the hot and cold temperatures are respectively set to 
be 893 K and 293 K. In order to using the same size stack and sub¬ 
sequently the same boundary layer condition for both working 
gases, Helium and Air, the frequency for Helium is chosen to be 
50 Hz which corresponds to 6.5 Hz for Air. The drive ratio, DR is 
supposed to be 5%. In fact, when the drive ratio is fixed, this leads 
to say that the size of heat exchangers is modified to match the de¬ 
sire drive ratio after finding the optimized values. 

The variable parameters are: 

The center of the stack x*, which varies between y and y + 0.5. 
By using dimensionless number, the wavelength is excluded from 
use in the calculation, which simplifies the comparison between 
both gases used in this paper. 

The hydraulic radius r h has three different values 0.15 mm, 
0.2 mm, and 0.31 mm, where the thermal penetration depth is in 
order of S t 0.3 mm. 

The traveling-standing wave ratio z has three different values 
0.6, 0.75 and 0.95. z depends on the boundaiy conditions of the 
resonator. So, after finding the optimized value of i, the boundary 
conditions will be chosen to match the desired z. 

After highlighting the data, Figs. 3 and 4 are resulting. 

As shown in Figs. 3 and 4, the acoustic power is strongly 
depending on the stack’s position, the stack’s hydraulic radius 
and the traveling-standing-wave ratio. As results show, at fixed z 
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Fig. 3. Synthetic optimization of rj ex vs. A(VV*) for Helium as working gas at 
frequency 50 Hz. 



Fig. 4. Synthetic optimization of rj ex vs. A(W*) for Air as working gas at frequency 
6.5 Hz. 


layer conditions. In other word, a wholesale machine with Air as 
working gas is 7.7 times longer than a machine working with He¬ 
lium. And we also note that the ratio of the acoustic power with 
Helium to the acoustic power with Air is equal to: 

/pgCo\ jf PgCo\ =15116 x10 9/ 6 1 2 70 x 10 8 _ 2 5 (27 ) 

V 7 / helium! V 7 ) air 

In other word, the acoustic power produced by the same stack 
using helium is 2.5 better than using the air. 


6. Optimization using the Particle Swarm method (PSO) 

The Particle Swarm Optimization (PSO) is an iterative method 
that tries to optimize a function or a problem. Every problem has 
multiple candidate solutions, called particles, which are character¬ 
ized by their positions and velocities over a mathematical formula, 
in a known search-space of dimension N. Every particle saves its 
own best position overall iterations and saves also the best posi¬ 
tion found so far in the search-space by which all the swarm (the 
group of particles) is expected to move toward it. 

The Particle Swarm Optimization (PSO) method was invented 
by Kennedy and Eberhart [28]. The PSO is then used when classical 
optimization methods such as gradient descend and quasi-newton 
cannot be applied because of the complexity, irregularity or non¬ 
linearity of the problem. In addition, the PSO needs just few or 
no assumptions to optimize any problem and it can be employed 
in a large search-space dimension problem. Because the PSO meth¬ 
od is a stochastic method, it may not converge to the desired opti¬ 
mal solution from the first simulation, and hence, it is better to 
repeat the optimization procedure more than twice under the 
same conditions to ensure that the function is converging to the 
right value. 


and r h , the acoustic power varies in an interval range having a max¬ 
imum and a minimum, depending on the stack’s position in the 
resonator. This explains the variation of the oscillating pressure 
p*(x*) and oscillating velocity u*(x*) inside the resonator. The 
length of the interval range of the acoustic power is changed with 
the variation of r and r h . For the highest acoustic power, as Figs. 3 
and 4 show, the traveling-standing-wave ratio must be near a pure 
traveling-wave when t is equal to 1, and the stack’s hydraulic ra¬ 
dius, r h , must be in the order of the thermal penetration depth 
r h a 5 t which is approximately equal to 0.3 mm. The exergetic effi¬ 
ciency that corresponds to the highest acoustic power, which is 
approximately equal to 0.0011, is quite low and it is in the order 
of r] ex k. 1%. 

The highest exergetic efficiency is obtained at a traveling- 
standing-wave ratio near to a pure standing-wave, that is for 
t = 0.5, and for a stack’s hydraulic radius less than the thermal pen¬ 
etration depth, i.e. for r h <c 5 t . In this case, the exergetic efficiency 
has also an interval range with maximum and minimum values 
that depend on the stack’s position in the resonator. But the length 
of this interval is tending to be zero, or in other words, the exerget¬ 
ic efficiency becomes independent of the stack’s position in the res¬ 
onator when r is a pure traveling-wave, i.e. when t = 1, and when 
r h is in the order of thermal penetration depth, i.e. r h r* c5 t . At fixed r 
and for r h varied increasingly, or at fixed r h and for t varied increas¬ 
ingly, the highest exergetic efficiency is decreased. The acoustic 
power corresponding to the highest exergetic efficiency, which is 
approximately equal to 75% for Helium and 50% for Air, is quite 
low and it is in the order of AW* « 0.0001. 

As we said previously, the frequency for helium and air as work¬ 
ing gases must be equal to 50 Hz and 6.5 Hz respectively in order 
to obtain the same stack size and subsequently the same boundary 


6.1. Algorithm 

We define the cost function to be optimized (minimized or 
maximized) as: 

r n -> OS 

where x is the position vector of dimension N. Here, N represents 
the dimension of the search-space. 

Let P be the number of candidate solutions (particles) in the 
swarm, where each particle has a position vector of dimension N. 
So we define x,;i = 1 ,...,P as the position vector of a particle i 
and dj as its velocity. 

Let k be the number of iterations, so the algorithm for maximiz¬ 
ing a function can be summarized as: 

I. for each iteration j = 1. k, 

1. for each particle i = 1 ,,P, 

i. Evaluate the fitness value of the cost function, jj'(xj) at 
position xj. 

ii. If fj(x ;:) >fi est (x'best)- then /U(*U) and 

XU = xj, where f best (x' best ) is the best value retained by 
a particle i at position xU- 

ill. If /'(xj) >/U(xLt). then /L t (*Lt) =/*(*;) and 

where /U(*U) * s the global best value retained at 
position xU. 

iv. Update particle position xj +1 and velocity z>j +1 . 

v. If stopping conditions are satisfied go to step II. 

II. Report results and terminate. 
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6.2. Defining cost functions to be optimized 

We are interesting by maximize three functions: 

• the first one is the exergetic efficiency of the thermoacoustic 
engine t) ex , 

• the second one is the acoustic power produced by the stack of a 
thermoacoustic engine, A{W), 

• And the last function is the product between the acoustic power 
produced by the stack and the exergetic efficiency, ri ex x A(W). 
This function is useful specially to find a tradeoff between the 
best of ;; ex and A(W* *). Because the last section has shown that 
when ri ex is in its maximum, A(W) is in its minimum, and vice 
versa, so, it is better to design a thermoacoustic device based on 
this function. 

6.3. Defining search-space and parameters of the problem 

For the first time, the PSO is involved in the thermoacoustic re¬ 
search, so in this paper we concern ourselves by studying a 3D 
search-space dimension problem. 

The cost functions studied in this paper depend on three 
parameters: 

• The hydraulic radius of a parallel-plate stack, r h , 

• The center of the stack in the resonator, x*, and 

• The traveling-standing-wave ratio, t. 

These three parameters are studied to see their influence on the 
maximum cost functions mentioned above, i/ ex , A(W*) and 
fcxA <W). 

The search-space of particles is summarized in Table 2. For x*, 
the interval range corresponds to the stack’s position varying from 
zero to the half of a wavelength. The hydraulic radius, r h , is lying 
between 0.15 mm and 0.31 mm where the thermal penetration 
depth, t> t , is in order of 0.3 mm, the second value is in order of 5 t , 
and the last value is greater than S t . The traveling-standing-wave 
ratio, t, will be varied between 0.5, which corresponds to a pure 
standing-wave ratio, and 1, which corresponds to a pure travel¬ 
ing-wave ratio. In this paper, we concern ourselves in what is hap¬ 
pening inside the stack of a thermoacoustic engine. Then, the 
boundary conditions and the size of the heat exchangers are cho¬ 
sen to match the optimized values obtained. 

Table 3 summarizes data of fixed parameters concerning cost 
functions for each of the two working gases studied, Helium and 
Air. The mean pressure of the system is supposed to be 25 bar, 
the hot and cold temperature are equal to 893 K and 293 I< respec¬ 
tively, while the dimensionless length of the stack is set to be 
0.0025. The drive ratio is fixed also at 5%. We suppose that the drive 
ratio is fixed in order to say in other words that the size of the hot 
and cold heat exchanger can be changed to maintain the desire va¬ 
lue of drive ratio. Thermophysical properties of helium and air used 
in this paper are also found in Table 3. Concerning the frequency, 
the thermoacoustic engine is working at 50 Hz when Helium is used 
and at 6.5 Hz when Air is used. The difference in the frequency be¬ 
tween Helium and Air aims at maintaining the same stack size and 
subsequently the same boundary layer conditions. 


Table 2 

Search-space for optimizing. 


Parameter 

Symbol 

(unit) 

Helium 

Air 

Dimensionless number of center of 
the stack 

K 

IM+o-sj 

5;f+°'5] 

Hydraulic radius 

r h (mm) 

[0.15: 0.31] 

(0.15; 0.31] 

Traveling-standing-wave ratio 

T 

10. 5; 1] 

[0.5; 1] 


Table 3 

Parameters concerning functions to be optimized. 


Parameters 

Symbol (unit) 

Helium 

Air 

Mean pressure 

Pg (bar) 

25 

25 

Hot and cold temperature 

T„. T 0 (K) 

893, 293 

893, 293 

Dimensionless number of length of 
stack 

L* 

0.0025 

0.0025 

Frequency 

fr (Hz) 

50 

6.5 

Drive ratio 

DR 

5% 

5% 

Dynamic viscosity 

p (kgms- 1 ) 

0.000018 0.0000175 

Gas constant 

R (J Kg- 1 K 1 ) 

2080 

287 

Thermal conductivity 

kg 

(Witt 1 K" 1 ) 

0.1392 

0.02415 

Constant pressure specific heat 

C p a Kg- 1 K- 1 

) 5183 

1005 

Specific heat ratio 

7 

1.667 

1.4 


Table 4 

PSO parameters. 

PSO parameters 

Maximum number of K = 2,000,000 

iterations ( K) 

Number of particles P= 24 

(P) 

Stopping conditions If, for 250 iterations, the variation of global best, 
/^(xf ) is less than l-e25, then stop the 
simulation 


Regarding the particle swarm optimization method, we will re¬ 
fer to the particle swarm optimization toolbox developed by Birge 
[29]. Table 4 summarizes initialization parameters of PSO toolbox 
used for our optimizations case. 


6.4. Optimization results and discussion 

The particle swarm optimization method is involved for the first 
time in the thermoacoustic research in this paper. Hence, to be 
confident that this method works well, two verification methods 
are done: 

• As the PSO method is a stochastic method, 10 trials were done, 
under the same conditions, for each optimization function in 
order to ensure that the results were converging to the right 
solution and to avoid the stagnation problem of such optimiza¬ 
tion method. 

• In addition, optimization of the exergetic efficiency and acoustic 
power are studied in order to compare them with the results 
existing in the literature and therefore to guarantee that the 
powerful PSO method works well. 

Then, as a novel contribution, the PSO method is used to opti¬ 
mize a third function for the first time in the literature, i.e. the 
product )j ex x A(W*}, the product between the exergetic efficiency 
and the acoustic power. 


6.4.1. Exergetic efficiency optimization i/ ex (r h ,x*,r) 

Figs. 5 and 6 show the evolution of the exergetic efficiency of a 
thermoacoustic engine working with Air and Helium, respectively. 
As results shown in Table 5, the best values of the exergetic effi¬ 
ciency of a thermoacoustic engine are 50.14% for Air and 74.94% 
for Helium. It is well known that the thermoacoustic devices work¬ 
ing with Helium are better than Air, because the thermophysical 
characteristics of Helium, which is a noble gas. The best values of 
the exergetic efficiency correspond to: 
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Exergetic Efficiency Optimization of a Thermoacoustic 
Engine working with Air at 6.5 Hz and DR = 5% 



- TRIAL 1 

- TRIAL 2 

TRIAL 3 

- TRIAL 4 

TRIAL 5 

- TRIAL 6 

TRIAL 7 
TRIAL 8 
TRIAL 9 
TRIAL 10 


Fig. 5. Evolution of the exergetic efficiency of a thermoacoustic engine working 
with Air. Here, */ ex (x* = 0.011;r* = 0.150 mm;i = 0.5) = 50.14%. 


Table 5 

Exergetic efficiency optimization results. 



Parameters 


Exergetic efficiency 

Acoustic power 


K 

r„ (mm) 

T 

lex (%) 

a(w; c ) 

For Air 

0.011 

0.0150 

0.5 

50.14 

0.0000824 

For Helium 

0.464 

0.150 

0.5 

74.94 

0.0000698 


0.001188 for Air and 0.001228 for Helium, which, by using Eq. 

A( vv* ) 

(27), are equivalent to A ^“" > « 2.5 \ _> A(W Hdium ) ss 2.6 

j 

xA (W/ar). These results are shown in Figs. 7 and 8 and are summa¬ 
rized in Table 6. 

The best values of the acoustic power correspond to: 


Exergetic Efficiency Optimization of a Thermoacoustic 
Engine working with Helium at 50 Hz and DR = 5 % 



Iteration 


Fig. 6. Evolution of the exergetic efficiency of a thermoacoustic engine working 
with Helium. Here, i7 ex (x* = 0.464:0, = 0.150mm;t = 0.5) = 74.94%. 


• Pure standing-wave, x = 0.5, which means, for example, a reso¬ 
nator tube with two closed end boundary conditions. 

• Hydraulic radius less than the thermal penetration depth, S h . In 
our case, r h = 0.5 mm, where <5 rs 0.3 mm. 

• Relative position of the stack in the resonator, x*, which is equal 
to 0.011 when Air is the working gas and 0.46 for Helium. 


The dimensionless acoustic power values corresponding to the 
best value of the exergetic efficiency are relatively low, 
82.4 x 1CT 6 for Air and 69.8 x 10~ 6 for Helium. If we are taking 
into account Eq. (27), we find that the ratio of the acoustic power 
for Helium to Air corresponding to the best value of exergetic effi¬ 
ciency is equal to 


A (VV Helium) 
A(W /lir ) 


2.5 


^(y^xc,Helium) 


a(W jMtam) « 2.1 X A (W Mr ) 


The obtained results, for the exergetic efficiency optimization, 
are logic and in accordance with those that can be found in the lit¬ 
erature. Kang et al. [23] shows that the best value of exergetic effi¬ 
ciency decreases as the traveling-standing-wave ratio increases as 
shown in this paper. Kang defines the pressure field in a way that 
the traveling-standing-wave ratio, t, is equal to 0 for a pure stand¬ 
ing-wave and 1 for a pure traveling-wave. Kang’s definition of x is 
different from that discussed in this paper, but the consequences 
on the exergetic efficiency are the same regardless of the represen¬ 
tation of the pressure field. 


6.4.2. Acoustic power optimization A(W* c )(r h ,x*, t) 

For the highest acoustic power, the dimensionless number of 
acoustic power’s best value obtained in this paper are equal to 


• Pure traveling-wave, t = 1, which means, for example, a circled 
resonator tube. 

• Hydraulic radius is in order of the thermal penetration depth, 5 t , 
in our case, r h rs 5 t rs 0.3 mm. 


Acoustic Power Produced Optimization of a 
Thermoacoustic Engine working with Air at 6.5 Hz and 


0.0012 


DR = 5% 



100 150 200 

Iteration 


250 


- TRIAL 1 

- TRIAL 2 

- TRIAL 3 

- TRIAL 4 

TRIAL 5 

- TRIAL 6 

TRIAL 7 

- TRIAL 8 

TRIAL 9 
TRIAL 10 


Fig. 7. Evolution of the acoustic power of a thermoacoustic engine working with 
Air. Here, A(W()(x* e = 0.2;r„ = 0.291 mm;t = 1) = 0.001188. 


Acoustic Power Produced Optimization of a 
Thermoacoustic Engine working with Helium at 50 Hz 
and DR = 5% 



Fig. 8. Evolution of the exergetic efficiency of a thermoacoustic engine working 
with Helium. Here, A(W()(x* E = 0.15;r h = 0.293 mm;r = 1) = 0.001228. 


Table 6 

Acoustic power optimization results. 



Parameters 

Exergetic efficiency 

Acoustic power 


*c 

r h (mm) 

X t] ex (%) 

A (W'xc) 

For Air 

0.2 

0.291 

1 1.07 

0.001188 

For Helium 

0.15 

0.293 

1 1.8 

0.001228 
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Table 7 

Acoustic power times exergetic efficiency optimization results. 



Parameters 

K 

r„ (mm) 

T 

Exergetic efficiency 

Oex TO 

Acoustic power 

a(vv; c ) 

Exergetic efficiency times Acoustic power 

lex X A(W( C ) 

For Air 

0.0516 

0.150 

0.61 

29.1 

0.0002166 

0.000063 

For Helium 

0.0587 

0.150 

0.63 

39.61 

0.0002494 

0.0000988 


• Position of the stack in the resonator, xj, which is equal to 0.2 

when Air is working gas and 0.15 for Helium. 

As Table 6 shows, the exergetic efficiencies, for both Air and He¬ 
lium, that correspond to the best value of acoustic power are quite 
low. These results are in conformity with the work of Kang et al. 
[23]. Kang found that the best value of the normalized acoustic 
power gain is when the traveling-standing-wave ratio t = 0.68. 
But this value is not matching the best value of the acoustic power 
itself. To calculate the acoustic power gain in the Kang’s paper, we 
must multiply the normalized acoustic power gain by 
p 2 (1 + t 2 ). If this transfer condition is taken into account, we 

ATpgC-T 3 v ' 

will find that the best value of the acoustic power gain is when 
the wave is a pure traveling, t = 1. 

6.4.3. Exergetic efficiency times acoustic power optimization 

In the previous simulations, we saw that the best value of exer¬ 
getic efficiency has an acoustic power relatively low, while the best 
value of the acoustic power is when the exergetic efficiency is very 
low. So it might not be useful to design a thermoacoustic device 
based only on the best value of the exergetic efficiency or on the 
acoustic power. Hence, the idea is to optimize the product between 
the exergetic efficiency and acoustic power which seems to be a 
tradeoff solution between the best values of the acoustic power 
and exergetic efficiency. Optimization results of this function are 
shown in Table 7. The best value of this product is when the trav- 
eling-standing-wave ratio is approximately 0.62. So, it is necessary 
to adjust the boundary condition of the thermoacoustic engine to 
match this traveling-standing-wave ratio. Figs. 9 and 10 show 
the evolution of the best value of the product of acoustic power 
and exergetic efficiency. 

Backhaus and Swift [14] have presented an empirical study in 
which they improved the exergetic efficiency of a thermoacoustic 
engine by making a traveling-standing-wave ratio between a pure 
traveling and standing wave and by making the stack’s hydraulic 
radius shorter than the thermal penetration depth. But they did 
not give any explanation about the value of the traveling-stand- 
ing-wave ratio neither about the optimized function. The mathe¬ 
matical optimization study presented in this paper has led to 
results (see Table 6) that are in conformity with those of [14]. 
Although this proves that Backhaus and Swift were in the right 
path, our mathematical study is general and can be applied to a 
wide range of empirical systems. Furthermore, this shows the 
importance of the contributions of this paper: the use of PSO and 
the optimization of the product function described before. 

7. Conclusions 

In this work, an optimization method, the particle swarm opti¬ 
mization method (PSO), for the thermoacoustic acoustic power and 
exergetic efficiency is represented and analyzed. It takes into ac¬ 
count the effect of the stack’s hydraulic radius, the stack’s position 
in the resonator and the traveling-standing-wave ratio. 

As a result, the optimal solution of the acoustic power is obtained 
when the wave is a pure traveling and when the hydraulic radius is in 


Exergetic Efficiency X Acoustic Power Produced 
Optimization of a Thermoacoustic Engine working with 



TRIAL 1 
TRIAL 2 
TRIAL 3 
TRIAL 4 
TRIALS 
TRIAL 6 
TRIAL 7 
TRIAL 8 
TRIAL 9 
TRIAL 10 


Fig. 9. The best value of the exergetic efficiency times acoustic power of a 
thermoacoustic engine working with Air. Here, A(W() x OexiKe = 0.0516; r h = 
0.150 mm;r = 0.61) = 0.000063. 


Exergetic Efficiency X Acoustic Power Produced 
Optimization of a Thermoacoustic Engine working with 
Helium at 50 Hz and DR = 5 % 



- TRIAL 1 

- TRIAL 2 

- TRIAL 3 

- TRIAL 4 

TRIALS 

- TRIAL 6 

TRIAL 7 

- TRIAL 8 

TRIAL 9 
TRIAL 10 


Fig. 10. the best value of the exergetic efficiency times acoustic power of a 
thermoacoustic engine working with Helium. Here, A:'W)) x t/„(x* e = 0.0587; o, 
= 0.150 mm;t = 0.63). 


the order of the thermal penetration depth. This leads to a circled res¬ 
onator tube with a stack having a hydraulic radius in the order of the 
thermal penetration depth. When the wave is a pure traveling wave, 
this means that the oscillating pressure and velocity are in phase 
which explains why the acoustic power is maximal when the wave 
is a pure traveling-wave. But, in this case, the exergetic efficiency 
is quite low because the thermal and viscous dissipations are high. 

For the highest exergetic efficiency, the thermoacoustic engine 
must work in a pure standing-wave, while the hydraulic radius 
should be less than the penetration depth which means that the 
thermal dissipation is low and this explains the highest efficiency. 
The resonator, in this case, has two closed end in order to obtain 
the standing-wave motion. Since the thermoacoustic engine, in 
this case, works at standing-wave, this leads to say that the oscil¬ 
lating pressure and velocity have 90° phase difference which gives 
explanation for the lowest value of the acoustic power. 

Because it is not much useful to build a thermoacoustic engine 
based only on higher exergetic efficiency with lower acoustic power 




























80 


H. Chaitou, P. Nika/Energy Conversion and Management 55 (2012) 71-80 


or on higher acoustic power with lower exergetic efficiency, the 
product between acoustic power and exergetic efficiency is studied 
also. This function is in fact more realistic, because it leads to obtain 
a good exergetic efficiency and acoustic power. The corresponding 
traveling-standing-wave ratio is in the order of 0.6 while the 
hydraulic radius is shorter than the thermal penetration depth. 

After obtaining interesting results by using the PSO method, our 
future work will focus on studying more complicated thermoacou¬ 
stic devices by trying to take into account a large search-space 
dimension problem. 
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